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Abstract 

We perform a detailed analysis of thermal leptogenesis in the framework of seesaw 
models which approximately conserve lepton number. These models are known to 
allow for large Yukawa couplings and a low seesaw scale in agreement with neutrino 
mass constraints, and hence to lead to large lepton flavour violating rates that can be 
probed experimentally. Although large Yukawa couplings lead to (inverse) decay rates 
much larger than the Hubble expansion rate, we show that the leptogenesis washout 
induced is generically small if the mass splitting between the right-handed neutrinos 
is small enough. As a result, large lepton flavour violating rates are compatible 
with successful leptogenesis. We emphasize that this scenario does not require any 
particular flavour structure. A small splitting is natural and radiatively stable in this 
context because it is protected by the lepton number symmetry. 



1 Introduction 

The seesaw models provide an attractive and straightforward explanation for both the 
recently observed tiny neutrino masses and the baryon asymmetry of the universe pQ. 
However, it is a very difficult task to probe these models experimentally, especially the 
most studied "type-I" seesaw model [2j based on the existence of gauge singlet right- 
handed (RH) neutrinos. In the generic situation, where lepton number violation and lepton 
flavour violation are generated at the same scale, neutrino mass constraints require either 
an experimentally unreachable heavy scale for the mass of the RH neutrinos, or tiny Yukawa 
couplings (leading to suppressed RH neutrino production cross sections). In both cases the 
rate of rare lepton processes, such as fi — )■ ey, is expected to be highly suppressed. 

Without assuming any extra new physics beyond the right-handed neutrinosj^] there 
exists however a theoretically motivated class of seesaw models where tiny neutrino masses 
do not imply that the lepton flavour violation is suppressed: the seesaw models which 
approximately conserve the total lepton number L [S]-[I2]- These models are based on the 
fact that, even in the presence of both RH neutrino masses and neutrino Yukawa couplings, 
it is possible to conserve L, and thus to keep the left-handed neutrinos massless. Hence the 
lepton flavour violating scale can be disconnected from the lepton number violating one (see 
[T2] for a detailed discussion). As a result, the Yukawa couplings are not constrained to 
be small even if the RH neutrino mass scale is low. Large lepton-flavour-violating (LFV) 

^blancheQumd.edu, thambye@ulb.ac.be, fxjossemichaux@gmail.com 

2 In the following we assume a non-supersymmetric setup. In supersymmetric seesaw models it is known 
that large flavour violating rates can be induced in agreement with leptogenesis, see e.g. OS]. 
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-but L-conserving- processes can be induced, with rates that could be observed if the RH 
neutrino mass scale is sufficiently low, a few tens of TeV or less. Neutrino masses can 
be explained subsequently by introducing small L-violating perturbations in the Yukawa 
coupling matrix and/or the RH neutrino mass matrix. 

Successful baryogenesis via leptogenesis is known to be feasible at the TeV scale [131 E] 
through the resonant leptogenesis mechanism [TB*l [TB| [T5] , which, through the virtual N 
propagator of the one-loop self-energy diagram (TTJ [TU EE] , requires a quasi-degeneracy of 
at least two of the RH neutrinos. A nice feature of the approximately L-conserving models 
is that they precisely yield a quasi-degenerate spectrum, with a mass splitting proportional 
to the small L-violating entries. Successful leptogenesis is achievable in this case, at least 
for Yukawa couplings suppressed enough [20J. 

In this letter we investigate the following question: is successful leptogenesis compatible 
with observable rare lepton processes? In Ref. [21] it was shown to be feasible considering 
a particular flavour structure (suppressed r Yukawa couplings) and with an extra 5*0(3) 
symmetry between the 3 RH neutrinos. In the following we show that it is generically 
possible even with two RH neutrinos and without the need to assume any particular flavour 
structure. This is due to the fact that the large Yukawa couplings which can lead to 
observable rare lepton processes are precisely the ones which do not break lepton number, 
therefore not necessarily causing any washout of the lepton asymmetry. 

In Section 2 we introduce the model with approximate L conservation which is going to 
be at the center of our study. In Section 3 we show explicitly how the washout is suppressed 
due to subtle interference effects in the L-conserving limit. In section 4 we compute the CP 
asymmetry and, in Section 5, study the consequences of our finding for the observability of 
rare lepton processes. In Section 6 we conclude and emphasize the generality of our results. 



2 The model 

We consider for illustration the simple Type A model considered in [12J with only two RH 
neutrinos, iV lj2 (see also [20J ) . With the charge assignment L — 1, L — — 1 to ^ and N 2 
respectively, the Lagrangian can be split into a L-conserving part and a L-violating one: 

C L = \N i 0N i -\Yj! 1 $n Loi + ^MN 1 Nl + h.cA , (1) 
Cf = - N 2 ft i La + ^N^t + ^ 2 e ia N 2 N< + h.cj . (2) 

with (ft = (0 + 0°) T the ordinary Higgs doublet. Note that in this Lagrangian yU 1; \i 2 and 
M are real parameters. The phases of /ii and M have been rotated away. The Yukawa 
couplings Y and the RH neutrino mass M do not break lepton number, i.e. do not induce 
any L-violating dimension 5 operator for neutrino masses. But they do generate a dimension 
6 L-conserving operator, and hence rare lepton processes, which can be large. For example 
a fi — > ey rate of order the experimental upper bound Br(/i — > cy) ~ 10~ n can be obtained 
for M ~ 1 TeV if the Y couplings are of order few 10~ 2 . On the other hand, the couplings 
Y', as well as the diagonal mass entries fi± and fi 2 , do break lepton number and are therefore 
constrained to be small by the neutrino masses: Y' <§C Y and /ii, fi 2 M. This can be 
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seen by diagonalizing the neutral state mass matrix which results from Eqs. 



Y?v Y'Jv 
M„ = \ Y a v iii M 
Y' a v M /i 2 e ia 



(3) 



and which, at first order in /i 2 and Y', gives the following light neutrino mass matrix 
(with v = 174 GeV): 



Tfly = V 



Y' T —Y + Y T —Y' 
M M 



Y T — a 2 e ia —Y 
M M 



(4) 



For M ~ ITeV a neutrino mass of order 0.1 eV typically requires Y'Y ~ 10~ 12 and/or 
/x 2 Y 2 /M ~ 10 -12 . For /i2 = this model has also the nice and rather unique property to 
allow for the full reconstruction of the flavour structure of the model from the values of 
the light neutrino mass matrix entries [12]. For yU 2 ^ it is not the case but still the full 
flavour structure of the L-conserving rare lepton processes can be reconstructed, leading to 
definite predictions for the various rare lepton processes, up to an overall normalization. 

To consider leptogenesis in this framework, it is convenient to go to the basis where 
the RH neutrino mass matrix is diagonal with real and positive entries. The Lagrangian is 
then given by 



i;Y, jJS, I h ia N^U La + l -MiNiN° + h.c. 



C 



To first order in /x/M and Y', the mass eigenvalues are given by 



1,2; a = e,/i,r). (5) 



Mi 



(6) 



where we defined /ii + /i 2 e ia = /xe 1 ^, so that /i 2 = /i 2 + /i 2 + 2/ii/x 2 cos a. As for the Yukawa 
couplings, again to first order in fi/M and Y', we find 



hi c 



h 



2a 



V2 
1 

7? 



-i(*-A)/2 



-i(0+A)/2 



1 + 



1 - 



where 



A = sin a 



f4~f4 

A* 2 ~ A* 1 ! 
A*1^2 



j*Y - YL 



^Y a + Y' 



(7) 
(8) 

(9) 



3 Washout from AL = 2 scatterings 

The RH neutrinos being unstable particles, one cannot in principle define asymptotic free 
states for them. However, as shown in (TSJ [22] , their interaction properties (decay, inverse 
decay...) can be deduced from the AL = 2 scatterings mediated by on-shell heavy neutrinos 
(■a^ ^ £p<t>- The corresponding 7^1=2 a interaction rate enters in the Boltzmann equation 
for the lepton flavour a as 

,n/,A f Y Nj(z) ,\ n Y e a (z) to t / 1f y, 

^H{z)— = ^ Cfa ^_-lj 7 „ < -^_2 7A£ ^ iaI (10) 
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where s is the entropy density, z = Mx/T, H(z) is the Hubble expansion rate, Yx = nx/s 
is the number density of particle X normalized to the entropy density and ej Q is the CP 
asymmetry from the RH neutrino N t for the flavour a. Note that we consider a Boltzmann 
equation for each flavour, even though our results will not require to assume any particular 
flavour pattern. 

One can separate the AL = 2 scattering into an on-shell and off-shell contribution, 
7al=2q = 7al=2q + 7ai=2 «■ As shown in the appendix of |23j, in usual leptogenesis 
scenarios, the on-shell part of the AL = 2 scattering, which is the dominant contribution, 
yields the inverse decays of leptons into heavy neutrinos, 

tSU* = % = = ^-^MK&Mim (n) 

where fCi^z) are the modified Bessel functions. Note that Ylia^N- a = 7jv Inverse decays 
contribute to a depletion of the lepton asymmetry, parametrized by the washout parameter 
Ki a = TNi, a /H(T = Mjvi), where the Ni decay width is Tx,, a = \hi a \ 2 Mi/8 ir. With 
quasi-degenerate RH neutrinos the relevant washout parameter is given by the sum K a = 
K Xa + K 2a = \Y a \ 2 M/8ir/H(M) + 0(Y' 2 ^ 2 /M 2 ). 

There is an apparent contradiction here. On the one hand, it is clear that in the 
L-conserving limit AL = 2 scatterings must vanish. On the other hand, the on-shell 
contribution quoted above clearly does not vanish, since the decay width is dominated by 
the L-conserving Y couplings, \h\ 2 ~ \Y\ 2 . For example, if the Y couplings are large, say 
10~ 2 , and if the RH neutrino mass M is as low as 1 TeV, the washout parameter K a is huge, 
about 10 9 , which gives a baryon asymmetry way too small. However, this huge washout is 
independent of the L-violating parameters, i.e. it would be present even in the absence of 
L-violation, which is clearly erroneous. 

In the following we show that this apparent contradiction disappears when one takes 
into account, on top of the inverse decay terms, the usually neglected interference terms 
involving different iVj's. To this end we start from the AL = 2 scattering reaction density 
as given from the reduced cross section d-& L=2ta : 



1 M 4 f°° 

7A°i=2,a(z) = --TTT-i \ dXy/xd& L=2 , a {x)K,i(zy/x) 
Z 04 71 J Xthv 



(12) 



where x = s/M^ , and with the scattering threshold given by x t hr = (mh + mi) 2 /M 2 ~ 0. 
The s-channel reduced cross section is given by^ 



&AL=2,a = \hla\ 2 A u ^ + l^2a| 2 ^22 ^ \hjf\ 2 + 2 Re 







(h la h* a ) A 12 ( h W h h) 
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,(13) 



where the h ai 's are the neutrino Yukawa couplings of Eqs. (J^J) — (J8j) - In the above expression 
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with Pr 1 = ( x — di + iy/EiCi) the dressed inverse propagators, whose poles are located 
at hi = cii — iy/afii, with a* = (Mj/Mi) 2 and q = (Ti/Mi) 2 . In the following, whenever the 
small L-violating c 2 — c% difference is irrelevant we will take c\ = c 2 = c. 



3 Here for simplicity we only consider the on-shell part, but it can be easily shown that the total AL = 2 
rate vanishes in the L-conserving limit, and therefore so does the off-shell part. 
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Using a narrow width approximation^] the AL = 2 washout term is given by 

Z 64 7T 4 47T ' 



7al=2 



X 





,2 i „2 



a{ + a 2 2 v /a 1 a 2 {y/aic 1 + yja 2 c 2 



y/c (a 2 - ai) 2 + (-y/olci + v/a 2 c 2 ; 



2 



f + af) + o(r 2 , / , 2 ) 



The first term in the parenthesis of Eq. (16) arises from the pure N\ and N 2 contributions, 



and yields inverse decay terms as in Eq. (11). The second term is an interference effect 



which is usually neglected in standard seesaw models. However, in the models we consider, 
this cannot be done, since both contributions are of the same order. 

Factoring out the L-conserving contribution to the 7^ inverse decay terms, which 
up to 0(Y' 2 , /z 2 /M 2 ) terms is given by a = (7^ a + jj} Q )/2, one finally obtains for the 
total washout 

on IN, a n (. . r r- 1 + 35^ x /^/2 „2\\ ^N,a 2 8 2 

= - r -2\l + 25Vi- 1 + SV - C + S2 +0<y ,»))= — ■ 1 + v - cS + 62 , (16) 

where the mass splitting enters in the parameter 8 = (M 2 — M\)jY\ = ///(M-y/c) ~ 
(a 2 — ai)/2\fc. There is an obvious cancellation in the parenthesis between the first two 
terms (pure N\ and N 2 contributions) and the third one (Ni-N 2 interference), leaving a 
washout which vanishes in the L-conserving limit, 8 — > 0. Therefore, no matter how huge 



the "naive" inverse decay term in Eq. (11) is, there is no washout in the L-conserving limit, 
as expected. Conversely, for a large mass splitting 5 ^> 1, the interference term becomes 
negligible, and one recovers the usual washout from inverse decays. This shows that the 
interference term which is rightly neglected in typical seesaw models (see e.g. [231 123]) can 
be essential in the context of approximate L-conserving frameworks. 

Practically, since the washout vanishes in the L-conserving limit, it remains small as 
long as the L-violating perturbations, Y', /i lj2 , remain small enough: it is suppressed as 
long as 5 <C 1. The relevant washout parameter including the interference term turns out 
to be 

X eS — K - 5 ~ l K ■ 5 2 (17) 

1 + y/cb + b 2 

For the numerical example above with M ~ 1 TeV and K a ~ 10 9 , if 8 < 10~ 3 one ends 
up with an effective washout ~ K a ■ 8 2 < 10 3 , which is suppressed enough to keep a 
substantial lepton asymmetry. 

To sum up, in the class of models we are considering, the washout of the asymmetry 
can be suppressed by many orders of magnitude, granted that the degeneracy parameter 
8 is small enough. A splitting <5 <§C 1 is technically natural since it is controlled by L- 
violating parameters. As a matter of fact, the radiative and thermal corrections to 8 are 
proportional to YY' and not to Y 2 as in typical seesaw scenarios. It is interesting to notice 
that, contrary to usual seesaw models with two RH neutrinos where the total washout is 



lower-bounded by the solar scale K\ + K 2 > K so \ ~ 9, the effective washout in Eq. (17) is 
lower-bounded only by terms quadratic in the small L-violating parameters Y' and fi/M. 
Therefore, Kf < LT so i can be achieved in this case. 



4 Or equivalently, up to negligible terms, taking the residue of the integral in Eq. (12) at the physical 
poles Si 2- 
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4 CP asymmetry 



The CP asymmetry generated during the decays of Ni into the lepton flavour a is given 
by 



3& y 



Im 


h* h ■ 

"'ia"'ja 




h* h ■ ^ 








2 



/: J + 



Im 


/?*/?■ ( 


^ '7 
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where /„ is the usual L- violating self-energy and vertex loop factor and f c is a L-conserving 
self-energy loop factor. In the limit we are interested in, namely M 2 — Mi, only the 
self-energy correction is relevant, and we get f 2,1 ~ — /*' 2 — f 2 ' 1 — —fc' 2 = /self with 
1221 



/self - " „ 2 ' ——2 ■ (19) 

{a 2 -a 1 ) + {^a 2 c 2 - y/aiCi) 

It is instructive to express the CP asymmetry in terms of the parameters of the original 
Lagrangian, Eqs. ([TJ-Q: 

\Y a \ 2 ( . w^^MW, .. 

ei« = e 2a ~ — sin at— — : H ^ — r— / se if . (20) 



4tt V 2fiM E/j'lVl 

The above expression shows explicitly that the CP asymmetry in each flavour crucially 
depends on the L-violating parameters, and therefore vanishes when L is conserved]^] 



At the singularity Mi = M 2 , the self-energy correction above, Eq. (19), has a reg- 
ulator given by y/a 2 c 2 — y/aiC\. In usual seesaw scenarios, one typically has \s/a 2 c 2 — 
sJaiC\\ ~ maxi (y'ajCj) because the decay widths r7v 12 are very different from each other. 
However, in our case, while the ^foiCi do not vanish in the L-conserving limit, y 'a 2 c 2 — y r a\Ci 
does. As a result \y/a 2 c 2 — yJaiC\\ <C min i ( v /ajCi). Therefore the resonance, which occurs 
for a 2 -ai = ^c 2 - ^c[, i.e. for 5 = 5 rcs ~ 2Re(FF , V < ^)/|F| 2 - (/i^-/if)/2M/i < 1, 
is obtained for a mass splitting much smaller than the decay width. Around the reso- 
nance, where the CP asymmetry reaches its maximum value e ~ 1, the CP asymmetry 
is no longer suppressed by the small L-violating parameters. In practice in the following 
we will never need to use such a small mass splitting, since in this case, unless the CP 
phases are highly suppressed, the baryon asymmetry produced by leptogenesis would be 
far bigger than the observed value. Therefore, we will always stick to the natural range 
\y/a 2 c 2 — sJa,iC\\ <C a 2 — a% <C max i ( v /a i Cj), or equivalently 5 res < i5< 1. In this regime, 
the self-energy contribution is simply given by 

* 2^' (21) 

It is remarkable that the conditions to get a large enough CP asymmetry and a suppressed 
washout are the same: 5<C 1. 



5 In models with more fermion singlets, a term in the self-energy diagram may remain in the limit of L 
conservation, see e.g. B — L conserving models of leptogenesis [26] or models with an extra N 3 [211 127] . 
Unless an extra symmetry beyond L conservation is assumed to have a quasi-degeneracy between at least 
3 fermion singlets, such a term turns out to be generically irrelevant for the generation of an asymmetry 
at the TeV scale due to lepton flavour equilibration [2"T] . 
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Note that, using a different method, a regulator oc s/OiCi was obtained in [HI EI], giving 
a resonance peak for 8 ~ 1. In the following we will stick to the regulator found in [22] 123]. 



Eq. (19), for two reasons. First, in [25], the self-energy contribution was calculated using 
two quantum field theory methods, and both converge to the regulator in Eq. (19). Second, 
the very different oscillation formalism typically applying to the K°-K and B°-B systems 
was used in PS] to compute the leptogenesis self-energy diagram in the limit of small mass 
splitting and the same result was obtained. Note also that the results of [23], Eq. (19), are 
obtained in the limit where the off-diagonal corrections to the propagator in the N1—N2 
system are small compared to the diagonal ones. More precisely, in terms of the parameters 
of the model, the condition of Ref. [23] reads [(tfh) 12 + (tfh) n ]/(32ir 2 ) < (M 2 - M x )/M 1 . 
In our framework, since the off-diagonal corrections are proportional to the L-violating 
parameters, i.e. {h)h)\2 + (h^h) 2 i oc 0(Y'Y, /j,/M), this condition is clearly fulfilled. 



5 Leptogenesis vs. rare lepton processes 



It is well known that for RH neutrino masses of order the weak scale and Yukawa couplings 
of order fewlO -2 , rates for rare lepton processes such as /i — )■ close to the present 
experimental bound, Br(/i — » erf) < 1.2 x 10 -11 [28], can be obtained [29] .[30]. From 



the suppression of the washout in Eq. (17) and the form of the CP asymmetry shown in 



Eqs. (20) and (21), it is easy to convince oneself that successful leptogenesis is possible 
at the same time. Taking for example M ~ 1 TeV, Y ~ few 1(T 2 , Y' ~ Yfx/M ~ 1(T 10 
(which may lead easily to the correct neutrino masses pattern), one gets 5 ~ 10~ 5 and thus 
a CP asymmetry e ~ few 10~ 6 . In this case, even though K a ~ 10 11 , the effective washout 
in Eq. (17) is reduced by a factor 1/5 2 ~ 10 10 , so that K cG ~ 10 2 . Such values of the 



washout and of the CP asymmetry are in the right ballpark to explain the observed baryon 
asymmetry of the Universe. 

To study this result more quantitatively it is convenient to use the Casas-Ibarra 
parametrization of the Yukawa couplings [31], which, in the 2 heavy neutrino case and 
for a normal hierarchy in the light sector, reads 



h 



2a 



(v / ^2 cos z U* 2 + y/rn^ sin zU* 3 ) , 
{—\frn2 sin z U* 2 + v /m^ cos z U* 3 ) 



(22) 
(23) 



The inverted hierarchy is obtained simply by exchanging the indices 2 — > 1 and 3 — > 2. In 
this parametrization, z = z a + iz^ is a complex angle, and U is the PMNS matrix. We will 
adopt the parametrization 



/ C12C13 S12C13 si 3 e 1UL1 \ 

U = -S12C23 - Ci 2 S23Si 3 e 1<5cp C12C23 - Si 2 S23Si 3 e 1<5cp s 23 Ci 3 •diag(l,e" 17/2 ,l), 
\ S12 S23 — C12 C23 S13 e 1<5cp — C12 s 2 3 — S12 C23 S13 e 1<5cp C23C13 J 

(24) 

where = sin#jj, Qj = cos 8ij. Neglecting experimental errors, we will use 
sin 2 (2#i 2 ) = 0.87, sin 2 (2# 23 ) = 0.92 [28], sin 2 (^ 3 ) < 0.03 (la) [32], whereas for the 
neutrino mass-squared differences we will take Am 2 ol ~ 7.59 x 10~ 5 eV 2 and Am 2 tm ~ 
2.43 x 10~ 3 eV 2 [28]. Recall that with only two RH neutrinos the lightest active neutrino 
is massless. Therefore, if the neutrino mass hierarchy is normal, one has mi = 0, 1712 = 



7 



a/ Amg 0l and m 3 = y/ Am s 2 ol + Am^ tm , whereas if it is inverted, m 3 = 0, m 2 = ■>/ Am^ tm 
and mi = v 7 Am ltm ~ 

In order to use the parametrization of Eq. ( |22| in a regime corresponding to the model 
we consider, Eqs. Q-(|2|, let us invert Eq. to write the couplings Y and Y' as 

/ M p Zb 

Y a = ^J--e~^{-i^- 2 U: 2 + ^U: z ) , (25) 



n = eW Vf ^ e " a + v^^s) + +iA) • (26) 



It is clear that the parametrization above yields Y 7 <C K if e -26 <C 1, the fi/M term being 
anyway suppressed by the approximate lepton number symmetry. In this limit, the CP 
asymmetry, Eq. ( 20 ) , takes the simple form: 

e a = 4 P a m3 ~ m2 sin (2z a ) ^ (27) 

m 3 + VTL 2 

where 

p _ m 3 \U a3 \ 2 + 2^m 3 m 2 \vL\{U* a2 U a3 ) + m 2 \U a2 \ 2 

m 3 + m 2 

is the branching ratio of the decay of Ni into the flavour a, with Y2 a = 1- violation 
in Eq. (27) is controlled by the "leptogenesis phase" z a and is maximal for z a = ir/A. It is 
interesting to notice that the CP asymmetry in the case of inverted hierarchy is suppressed 
compared to normal case by a factor a/ Am s 2 ol / Am^ tm . As for the K^f washout parameter, 
Eq. (17), for 5 -C 1 it is simply given by 

K es „ M \Y a \ 2 2 m 3( 2)+m 2(1) 22 2 ^ 2 2 

where ~ 1.08 x 10~ 3 eV and Xi — 27 (45) for normal (inverted) hierarchy. Had we 
neglected the interference term, we would have obtained a huge washout, K a ~ xi e 2zb - 
Though the suppression of the washout oc 5 2 is crucial, it turns out that we will always 
remain in the regime of "strong" washout, defined as ^3. In this regime, where 
the dependence on the initial conditions disappears, and where subtle effects like thermal 
corrections are irrelevant [23], we obtain finally for the baryon asymmetry 

lO f \ p-4.32 2,, 

Y b = 37 E °- 5 7j^ » M ) * E P «° 16 -jiSr sin ( 2 ^) > (3°) 

en \ ex / \ a / 

where X2 — 8.7 x 10~ 5 (5.4 x 10~ 7 ) for normal (inverted) hierarchy. There is a mild 
dependence on the parameters of the PMNS matrix in -^a -0 ' 16 ' which, barring large 
cancellations, varies between 3.7 and 4.8 for normal hierarchy and 3.6 and 4 for inverted. 
In the following we will constrain the baryon asymmetry to lie within the 3a range of 
WMAP5, i.e. 8.1 x 10~ n <Y B < 9.5 x 10~ n [33]. 

Let us now turn to the low-energy observables. The branching ratio for the rare lepton 
process /i — > cy is given by |29j: 

huh*^ h 2e h 2fl 2 



Br(/i->e 7 ) ~ ^( 2v2 ) 2 



+ 



2 

i 1 xvi. 

3a ( m 3 + m 2 



32ti \ M 



M{ M 2 

2 \9e,\ 2 e^, (31) 
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where a is the fine structure constant, and 



9efi — ■ 

m 3 + m 2 

As before, the replacement 3 — > 2 and 2 — )■ 1 should be made for the case of inverted 
hierarchy. The dependence on the parameters of the PMNS matrix is clearly larger for 
the branching ratio than for the baryon asymmetry. Here, \g e ^\ 2 varies from 0.01 to 0.2 
for the normal case, and from 0.06 to 0.7 in the inverted one, also granted that no large 
cancellations occur. 



Using Eq. (30) successful leptogenesis implies a relation between z 6 and the degeneracy 



parameter S, which we use to rewrite the branching ratio as 



Br(u^e 7 ) ~ ^( m S + m 2 V ., / \ 2 ! V„ /'„ n l,i ) sin (2.;,, A 

19 /250GeV\ 2 fl0- 4 \ 31 , NOq o 
^ X3 x 10~ 12 I - j x ( — ) xsin(22 a ) - 93 , (33) 

with X3 taking values between 4.1 x 10 _1 and 7.6 for the normal case and between 4.2 x 10~ 2 
and 2.9 x 10 _1 in the inverted one. From the above equations it is clear that the smaller 
the degeneracy parameter 8 is, the bigger the Yukawa couplings can be without spoiling 



the success of leptogenesis, hence the larger the LFV rate [i — \ is. Eq. (33) shows that 
for M ~ 250GeV-lTeV, observable fx — > rates typically require 5 ~ 10^-10 _5 . This 
corresponds to a mass splitting (M 2 — M\)/M\ ~ 10~ 8 -10 -10 . Such a level of degeneracy 
is comparable to the one needed in usual resonant leptogenesis scenarios. 



To check that our analytical result of Eq. (33) is correct, it is necessary to integrate 
numerically the set of Boltzmann equations, as given in Eq. ( 10 ). Fig. 1 shows our numerical 
results for the simple case fii — 0, M — 250 GeV, 5qp = 7 = and #i 3 = 10°, taking for 
the leptogenesis phase z a its maximum value 7r/4. In this simple case, as fi± = 0, the 



phase a is unphysical and only the second term of Eq. (20) contributes. The parameter fj,% 
is directly responsible for the mass splitting between the RH neutrinos, while Zb sets the 
scale of the Yukawa couplings. Fig. 1 displays, as a function of these 2 parameters, the 
iso-contours of the LFV rate compatible with successful leptogenesis for both normal and 
inverted neutrino mass hierarchies. A sphaleron freeze-out cut-off has been applied at the 
temperature T ~ 130 GeV taking the Higgs boson mass rrih = 120 GeV [35]. We find a 



good agreement with the estimate of Eq. (33), up to the sphaleron freezeout which has not 



been taken into account in this equation. Also given in Fig. 1 are the LFV contour plot as 
a function of the parameters K eS and 5. We observe that the branching ratio reached can 
saturate or even exceed the present experimental limit, By(/i — > e-y) < 1.2 x 10~ n , depicted 
as a green line in Fig. 1. 

The // — > cy decay is not the only flavour violating process whose current experimental 
upper bound could be saturated in agreement with leptogenesis constraints. This is also 
possible for the r radiative decays, which experimentally are constrained to the bounds 
Br(r ->■ /ry) < 4.5 • 10~ 8 and Br(r ->■ e>y) < 1.1 • 10~ 7 [28J. It turns out that the results 
of Fig. 1 can be applied to these rates for a degeneracy parameter 5 about one order of 
magnitude smaller than for the fi — > c-f rate^J As explained in Ref. p2] the observation 
of several radiative decays would allow to overconstrain this model. In the type-I seesaw 



6 Thc Yukawa couplings in the \i-e channel should be somewhat suppressed with respect to the ones of 
the other channel(s) for the more stringent /i — > ej constraint to be satisfied. 
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Log 10 B R (ji^ e 7) Log 10 B R (ji-f e y) 




Figure 1: Contours of Br(/i — > ej) compatible with successful leptogenesis plotted against 
/i 2 and z b (left panel), and K eS and 5 (right panel). The upper (lower) values correspond 
to the normal (inverted) hierarchical neutrino spectrum. The green line shows the present 
experimental limit Br(/i — > cy) < 1.2 x 10~ n . 



models, other channels like [l — > eee or r — > 3/ are expected to be more suppressed than 
the radiative decays, see e.g. Ref. [30J. 



6 Softly broken L case: = 

At the end of the previous section, by setting /i X to in Fig. 1, we numerically considered a 
simple case where the CP asymmetry is proportional to the Y' couplings. It is interesting 
to discuss what happens if we take Y' — instead. This is the situation of the usual inverse 
seesaw models [5]- [8] where L is assumed to be softly broken, and where in full generality 
n pairs of N\ } 2 are usually assumed. In such a framework the Yukawa coupling matrix Y 
in Eqs. ([if- ([3]) is an n x 3 matrix, /i lj2 are n x n matrices and Y' is a null nx3 matrix. 
It must be stressed that in this case too leptogenesis can be compatible with observable 
flavour violation rates. This can be seen from the n = 1 case above which, through the 



first term of Eq. (20), gives the CP asymmetry 



\Y a \ 2 . n <5<1 \Y a \ 2 ^1^2 

eia = £2a ^ — sin a —— / self ~ — — — sin a , (34) 

47r 2fj,M 167T ytr 



where Eq. (21) was used in the last equality. It is worth mentioning again that the CP 
asymmetry is not suppressed by the small L-violating entries in the limit of small splitting 
5^1. The larger are the Y a couplings, the larger is the /i — > e^y rate and the larger is the 
CP asymmetry. Note that a non-vanishing value of both fii and /i 2 is required in order to 
have a non-zero CP asymmetry, since otherwise the phase a can be rotated away. As for 



the washout, it is again suppressed as soon as 5 <C 1, Eq. (16) 



The generalization to the inverse seesaw case with n > 1 pairs of N12 is straightforward. 
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In this case the N±s and A^'s form n pairs of quasi-degenerate states and, as in Eqs. (16) 
and (34), each pair gives an unsuppressed CP asymmetry and a small washout as soon 
as the mass splitting is smaller than the decay width. Note that when Y' a = 0, n > 2 is 
required in order to give at least 2 neutrino masses. To consider this case in full details is 
beyond the scope of this paper because it involves many more parameters. 



7 Conclusion 

We have shown that, in the context of approximately L-conserving seesaw models, it is 
rather easy to generate large lepton flavour effects in agreement with successful leptogenesis. 
These models provide all the necessary ingredients. They involve large Yukawa couplings 
without leading to large neutrino masses, as required to induce large rates. They predict 
a small right-handed neutrino mass splitting protected from large radiative corrections by 
the approximate L-symmetry. If the splitting is sufficiently small, 5 <C 1, the washout is 
suppressed, proportionally to the small L-violating entries. And, finally, from the same 
condition, 5 <C 1, the suppression of the CP asymmetry by the small L-violating entries 
is compensated by its resonant behaviour. In this way the experimental upper bound of 
any of the rare lepton processes /i — > ej, r — > /i7, r — > can be saturated in agreement 
with successful leptogenesis. We emphasize that this pattern can work with only 2 RH 
neutrinos and that moreover it does not require to assume any flavour symmetry. Only 
a lepton number assignment for the various RH neutrinos has to be assumed for lepton 
number to be approximately conserved. 

To illustrate these results we have considered a minimal 2 RH neutrino model which 
has the additional virtue that from the knowledge of the neutrino mass matrix one can 
predict the rates of any rare lepton process up to an overall normalization scale. Therefore 
this model can be overconstrained. To our knowledge, if one does not assume any extra 
particles beside the RH neutrinos, there is no neutrino mass and leptogenesis seesaw model 
more testable than this one. 

Our results also apply to the usual inverse seesaw models where L is softly broken, with 
n pairs of fermion singlets and Y' = 0. Similarly it also applies to the 3 RH neutrino case 
where, if L conservation is approximate, two RH neutrinos are naturally quasi-degenerate. 
In particular they apply to the 3 RH neutrino models considered in Refs. [91 [TUl 112] . For 
the third RH neutrino, one has just to make sure that its mass is either sizably heavier 
than the two others or that at least one of its Yukawa coupling is suppressed enough. 

Finally it is worth mentioning that our main result about the washout suppression 
also applies to the Type-Ill seesaw models with the same Yukawa coupling matrix and 
heavy state mass matrix. However, leptogenesis suffers in this case from the high degree of 
thermalization of the fermion triplets by the gauge interactions [36] . It would be interesting 
to study more extensively this case to determine which region of the parameter space is 
allowed. 
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